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AERONAUTIC SYMBOLS 


1. FUNDAMENTAL AND DERIVED UNITS 


-\ 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbreviation 

T iPn rrfh 

l 

meter 

m 

foot (or mile) 

ft (or mi) 

JLJV/JLig l/Li — — — — — — 

Timp 

t 

second 

s 

second (or hour) 

sec (or hr) 

Force 

F 

weight of 1 kilogram. _ _ _ 

kg 

weight of 1 pound 

lb 

Power 

P 

horsepower (metric) _ 


horsepower _ 

hp 



/kilometers per hour 

kph 

miles per hour _ _ ... 

mph 

Speed-- 

V 

\meters per second 

mps 

feet per second. _ 

fps 


2. GENERAL SYMBOLS 


w 

9 

m 

I 


S 

S„ 

G 

b 

c 

A 

V 

2 

L 

D 

D 0 

D t 

D b 

c 




Weight =mg v 

Standard acceleration of gravity=9.80665 m/s 2 p 


Kinematic viscosity 
Density (mass per unit volume) 


or 32.1740 ft/sec 2 

W 


Mass= 

9 

Moment of inertia =mfc 2 . (Indicate axis 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


of 


Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft -4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, g 
True air speed 
Dynamic pressure, \pV 2 

Lift, absolute coefficient 

Drag, absolute coefficient 


3. AERODYNAMIC SYMBOLS 

i w Angle of setting of wings (relative to thrust line) 

i t Angle of stabilizer setting (relative to thrust 

line) 

Q Resultant moment 

12 Resultant angular velocity 


D 0 


Profile drag, absolute coefficient Dd 0 = ^ 


Di 


Induced drag, absolute coefficient C D .—^g 


Dry 


Parasite drag, absolute coefficient C Dp = ^ 
Cross-wind force, absolute coefficient 


C 


R 


a 

€ 

«0 

Oti 

<*a 


Reynolds number, p ^ where l is a linear dimen- 
sion (c.g., for an airfoil of 1.0 ft chord, 100 
mph, standard pressure at 15° C, the corre- 
sponding Reynolds number is 935,400; or for 
an airfoil of 1.0 m chord, 100 mps, the corre- 
sponding Reynolds number is 6,865,000) 
Angle of attack 
Angle of downwash 
Angle of attack, infinite aspect ratio 
Angle of attack, induced 

Angle of attack, absolute (measure J from zero- 
lift position) 

Flight-path angle 


Date Due 


It 26 ’5 

* 


t • 

l n 191 

& 



196 

) 



jpy 

1^1980 
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A SIMPLIFIED METHOD OF ELASTIC-STABILITY ANALYSIS FOR THIN CYLINDRICAL SHELLS 

By S. B. Batdorf 


SUMMARY 

This paper develops a new method for determining the buckling 
stresses oj cylindrical shells under various loading conditions. 
For convenience of exposition , it is divided into two parts. 

In part I, the equation for the equilibrium oj cylindrical 
shells introduced by Donnell in NACA Report No. Jf.79 to jind 
the critical stresses oj cylinders in torsion is applied to jind 
critical stresses jor cylinders with simply supported edges under 
other loading conditions. It is shown that by this method solu- 
tions may be obtained very easily and the results in each case may 
be expressed in terms oj two nondimensional parameters , one 
dependent on the critical stress and the other essentially deter- 
mined by the geometry oj the cylinder. The injluence oj boundary 
conditions related to edge displacements in the shell median 
surjace is discussed. The accuracy oj the solutions jound is 
established by comparing them with previous theoretical solutions 
and with test results. The solutions to a number oj problems 
concerned with buckling oj cylinders with simply supported 
edges on the basis oj a unified viewpoint are presented in a 
convenient jorm jor practical use. 

In part II , a modijied form oj Donnell 1 s equation jor the 
equilibrium oj thin cylindrical shells is derived which is equiv- 
alent to Donnell 1 s equation but has certain advantages in physical 
interpretation and in ease oj solution , particularly in the case 
oj shells having clamped edges. The solution oj this modijied 
equation by means oj trigonometric series and its application to 
a number oj problems concerned with the shear buckling stresses 
oj cylindrical shells are discussed. The question oj implicit 
boundary conditions also is considered. 

INTRODUCTION 

The recent emphasis on aircraft designed for very high 
speed has resulted in a trend toward thicker skin and fewer 
stiffening elements. As a result of this trend, a larger fraction 
of the load is being carried by the skin and thus ability to 
predict accurately the behavior of the skin under load has 
become more important. Accordingly, it was considered 
desirable to provide the designer with more information on 
the buckling of curved sheet than has been available in the 
past. In carrying out a theoretical research program for this 
purpose, a method of analysis was developed which is be- 
lieved to be simpler to apply than those generally appearing 
in the literature. The specific problems solved as a part of 
this research program are treated in detail in other papers. 
The purpose of this paper, which is discussed in two parts, 
is to present the method of analysis that was developed to 
solve these problems. 


In part I, the stability of a stressed cylindrical shell is 
analyzed in terms of Donnell's equation, a partial differential 
equation for the radial displacement w, which takes into 
account the effects of the axial displacement u and the cir- 
cumferential displacement v. Part I shows the manner in 
which this equation can be used to obtain relatively easy 
solutions to a number of problems concerning the stability of 
cylindrical shells with simply supported edges. The results 
of the solution of this equation are shown to take on a simple 
form by the use of the parameter k (similar to the buckling- 
stress coefficients for flat plates) to represent the state of 
stress in the shell and the parameter Z to represent the 
dimensions of the shell, where Z is defined by the following 
equations: 

For a cylinder of length L 



and for a curved panel of width b 


where 

r radius of curvature 
t thickness of shell 
and 

/x Poisson’s ratio for material 

The accuracy of Donnell’s equation is established by compari- 
sons of the results found by its use with the results found by 
other methods and by experiment. 

In the simplest method that has been found for solving 
Donnell’s equation, the radial displacement w is represented 
by a trigonometric series expansion. This method can be used 
to great advantage for cylinders or curved panels with simply 
supported edges but leads to incorrect results when applied 
uncritically to cylinders or panels with clamped edges. 

In part II an equation is derived which is equivalent to 
Donnell’s equation but is adapted to solution for clamped as 
well as simply supported edges by means of trigonometric 
series. This modified equation retains the advantages of 
Donnell’s equation in ease of solution and simplicity of re- 
sults. The solution of the modified equation by means of the 
Galerkin method is explained, and the results obtained by 
this approach in a number of problems concerned with the 
shear buckling stresses of cylindrical shells are given in 
graphical form and discussed briefly. Boundary conditions 
implied by the method of solution of the modified equation 
are also discussed. 
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SYMBOLS 


a length of curved panel (longer dimension) 

b width of curved panel (shorter dimension) 

d diameter of cylinder 


n, p, q\ 


V 

r 

t 

u 


w 


X 

V 

(L/titi) b mn , 
Cmri) d m n 


integers 

lateral pressure, positive inward 
radius of cylindrical shell 
thickness of cylindrical shell 

displacement in axial (x-) direction of point on 
shell median surface 

displacement in circumferential ( y -) direction of 
point on shell median surface 
displacement in radial direction of point on shell 
median surface; positive outward 
axial coordinate 
circumferential coordinate 


numeri cal co e ffici en ts 


shear-stress coefficient ^ T jj^> for cylinder or T jj\ for 
curved panel or infinitely long curved strip^ 
axial compressive-stress coefficient ( ^f)^ 


for 


c T x tb 2 


C 

1 V 

W {) 

D 

E 

F 


cylinder or for curved panel or infinitely 
long curved strip ^ 

circumferential compressive -stress coefficient 
i for cylinder or for curved panel or in- 
finitely long curved strip^ 

hydrostatic-pressure coefficient 

amplitude of deflection function 

plate flexural stiffness per unit length J 

Young’s modulus 

Airy’s stress function for the median-surface 

stresses produced by the buckle deformation 

/c ) 2 F . . • i r VF . . . 

, stress m axial direction; stress m cir- 


d 2 F 


cumferential direction; — shear stress 

L length of cylinder 

Q,Qi,Qy mathematical operators 


) 


Z 


curvature pan 
b- 


( J2 

^ Vl— ju 2 for cylinder or 


0 

X 


' A l — F 2 for curved panel or infinitely long 
curved strip^ 

L/\ for cylinder or b/\ for infinitely long curved 
strip 

half wave length of buckles; measured circumfer- 
entially in cylinders and axially in infinitely long 
curved strips 


£ dimensionless axial coordinate (x/b) 

t? dimensionless circumferential coordinate (y/b) 

H Poisson’s ratio 

t applied shear stress 

r CT critical shear stress 

g x applied axial stress, positive for compression 

g v applied circumferential stress, positive for compres- 

sion 

R s shear-stress ratio; ratio of shear stress present to 
critical shear stress when no other stress is acting 
U x axial-compressive-stress ratio ; ratio of direct axial 
stress present to critical compressive stress when 
no other stress is acting 


V 4 

operator ^ 

/y, &L Y_*. , o _*_+ *V 

\&e 2 ' by 2 ) dx 4 ^ dx 2 dy' 2 ' by* / 

X7 4 

v a 

operator ^ 


V 8 

operator ^ 

:(£+&)*) 

V G 8 

operator ^ 

:($+&)•) 

V" 4 

inverse operator defined by equation 


(V“ 4 (V 4 /) = V 4 (V" 4 /) =/) 

1. DONNELL’S EQUATION 

THEORETICAL BACKGROUND 

In most theoretical treatments of the buckling of cylin- 
drical shells (see references 1 to 3) three simultaneous partial 
differential equations have been used to express the relation- 
ship between the components of shell median-surface dis- 
placement u , v, and w in the axial, circumferential, and 
radial directions, respectively. No general agreement has 
been reached, however, on just what these equations should 
be. In 1934 Donnell (reference 4) pointed out that the 
differences in the various sets of equations arose from the 
inclusion or omission of a number of relatively unimportant 
terms (referred to in the present paper as higher-order 
terms), and proposed the use of simpler equations in which 
only the most essential terms (first-order terms) were 
retained. The omitted terms were shown to be small, and 
thus the simplified equations to be applicable, if the cylinders 
have thin walls and if the square of the number of circum- 
ferential waves is large compared with unity. Donnell 
further showed that the three simplified equations can be 
transformed into a single eighth-order partial differential 
equation in w (see appendix A of the present paper) in which 
the effects of the displacements u and v are properly taken 
into account; this equation will hereinafter be referred to as 
Donnell’s equation. 

When higher-order terms are included in the three partial 
differential equations previously mentioned, the resulting 
theoretical buckling stresses are usually very complicated 
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functions of the cylinder dimensions and the elastic proper- 
ties of the material. A family of curves is ordinarily drawn 
giving the critical stress as a function of the length-diameter 
ratio for specified values of the radius- thickness ratio and 
for given elastic properties (references 2, 3, and 5). When 
the higher-order terms are omitted from the equations and 
the requirements of an integral number of circumferential 
waves is removed, new parameters can be introduced which 
combine the cylinder dimensions and material properties in 
such a way that the results can be given in terms of a single 
curve. These parameters have been used, with slight varia- 
tions in detail, by Donnell, Kromm, Leggett, and Rcdshaw 
(references 4 and 6 to 9). The omission of the higher-order 
terms also greatly simplifies the calculations, and the calcula- 
tions are simplest if Donneirs equation, rather than the set 
of three simultaneous equations, is employed. DonnelPs 
equation, or an equivalent equation, may therefore be 
presumed to be the most promising for use in solving hitherto 
unsolved problems in the stability of cylindrical shells. 

In spite of the fact that it was introduced some time ago, 
DonnelPs equation has not achieved the wide acceptance 
for use in the stability analysis of cylindrical shells which 
it appears to merit. Some investigators have continued to 
use simultaneous differential equations in which higher- 
order terms appear, presumably on the assumption that the 
errors arising from neglect of these terms might lie undesir- 
ably large. Others have dropped second-order terms but have 


continued to employ simultaneous equations, probably in 
order to specify directly edge-restraint conditions having to 
do with displacements in the axial and circumferential 
directions, which cannot be done with DonnelPs equation. 

The purposes of part I are to establish the accuracy of the 
equation by comparing the results found by ' the use of 
DonnelPs equation with the results found by other methods 
and with experimental results and to investigate the question 
of boundary conditions on u and v. The additional purpose 
is achieved of presenting the solutions of a number of 
problems concerned with buckling of cylinders with simply 
supported edges on the basis of a unified viewpoint and in a 
convenient form for practical use. 

BUCKLING STRESSES OF CYLINDERS WITH SIMPLY 
SUPPORTED EDGES 

Lateral pressure. — The theory for the lateral pressure 
(uniform external pressure applied to walls only) at which 
a cylinder will buckle is given in appendix B in which it is 
assumed that the lateral pressure causes the buckling by 
producing a circumferential stress a y and that it affects the 
buckling in no other way. The results are shown in a 
logarithmic plot in figure 1. The ordinate in this figure is 
the stress coefficient k y which appears in the flat-plate 
buckling equation (see, for example, reference 3, p. 339) 
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(The discussion given in the section of the present paper 
entitled “Parameters Appearing in Buckling Curves” shows 
the relationship between a cylinder of length L and an 
infinitely long flat plate of width b=L.) The abscissa 



may be regarded either as a measure of the curvature, or, 
for any given ratio of radius to thickness, as a measure of the 
length-radius ratio of the cylinder. Figure 1 shows that 
for small curvature k v approaches the value 4, which applies 
in the case of simply supported long flat plates in longitudinal 
compression (reference 3, p. 327). As the curvature param- 
eter Z increases, the stress coefficient k y also increases. 
For large values of Z, the curve approaches a straight line of 
slope 1/2. This straight line is expressed by the formula 

ky= 1 .04Z 1/2 

As the length-radius ratio increases, for a given value of 
r/t, the number of circumferential waves n diminishes. Al- 
though n must be an integer, the curves of figure 1 were 
obtained on the assumption that n is free to vary continu- 
ously. Only small conservative errors are involved in this 
assumption. Because n= 1 corresponds merely to a lateral 
displacement of the entire circular cross section, the minimum 


value of n is 2, which corresponds to deformation of the 
section into an ellipse. This limitation on n results in 
splitting the curve of figure 1 into a number of curves for 
different values of rjt when Z becomes large. A cylinder 
r 

having a value of -=20 buckles into an ellipse when Ljr is 


about 10, and the value of Ljr at which such buckling occurs 
increases with increasing r/t. 

In figure 2 the curve of figure 1 is compared with results 
based on more complicated calculations given in reference 3 
and in reference 5. At fairly large values of Z the results 
given in reference 3 and in reference 5 are in good agreement 
with the results of the present paper. At small values of Z 
the curve based on reference 3 (Timoshenko) is definitely too 
low, because k v should approach the flat-plate value of 4 as Z 
approaches zero. An interesting feature of the comparison 
is that one calculation gives results below, and the other 
calculation results above, those given herein. The test data, 
taken from reference 5, are in reasonable agreement with 
and show more scatter than the theoretical curves. 

In the case of cylinders so long thatn=2, the requirement 
for the validity of Donnell’s equation that n 2 ^ 1 is no 
longer satisfied and appreciable error is to be expected. 
Indeed it may be shown that for very long cylinders when 
n— 2 Donnell’s equation gives 4 D/r 3 as the critical value of 
the applied lateral pressure, whereas the accepted theoretical 



Figure 2. — Comparison of present solution for critical circumferential-stress coefficients for simply supported cylinders with other theoretical solutions and with test results. 

(Timoshenko’s solution is from reference 3 and Sturm’s data and solution are from reference 5.) 
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result is SD/r 3 (by use of the formula given on p. 450 of 
reference 3). The curves for n = 2 will probably not often 

be needed, however, since they apply only wlicn^y^ 

which in the case of thin cylinders corresponds to a very 
large length-radius ratio, and if needed, the curves for 
n= 2 can be applied in conjunction with a correction factor 
0.75. 

Axial compression.— The theory for the axial stress at 
which a cylinder will buckle is given in appendix B, and the 
results are shown in figure 3. The ordinate is analogous to, 
and the abscissa identical with, the corresponding coordinates 
used in figure 1. Figure 3 shows that for small values of 
Z, k x approaches the value 1, which applies in the case of 
long flat plates in transverse compression with long edges 
simply supported (reference 3). For large values of Z, the 
curve becomes a straight line of slope 1. This straight line 
is expressed by the formula 

= 4 Z =0.7 0 2Z 

7T“ 


For any fixed value of rjt some value of Z always exists 
above which Ljr is so large that the cylinder fails as an 
Euler strut rather than by buckling of the cylinder walls. 
Pin-ended Euler buckling of cylinders is indicated in figure 3 
by means of dashed curves. 

The result just given for the critical-stress coefficient for a 
cylinder in axial compression leads to t he following expression 
for the critical stress: 


= 1_ 

V3~(W) 


Et 

V 


( 1 ) 



Figure 3.— Critical axial-stress coefficients for cylinders with simply supported edges. 
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The value given in equation (1) for the critical stress of a 
moderately long cylinder in axial compression by use of 
DonnelPs equation is identical with the value found by a 
number of investigators using other equations as starting 
points (references 1 to 3). In the case of cylinders under 
axial compression the errors involved in dropping the second- 
order terms are therefore concluded to be small. 

The buckling stresses given by equation (1) are neverthe- 
less in serious disagreement with the buckling stresses 
obtained by experiment (reference 10). For a discussion of 
the degree of correlation that can be found between theory 
and experiment for cylinders under axial compression, see 
reference 1 1 . 

Hydrostatic pressure on closed cylinders. — When closed 
cylinders are subjected to external pressure, both axial and 
circumferential stress are present. The theory for buckling 
under these combined loads is given in appendix B. The 
results are shown in figure 4. The ordinate C v used in this 
figure is a nondimensional measure of the pressure j) de- 
fined as follows: 


C v = 


prU 

7 t 2 D 


The coefficient C p can be directly related to the corresponding 
stress coefficients k x and k y . By definition 


, _(T v tL 2 
tt 2 D 


and, according to the hoop-stress formula, 



It follows from the three preceding equations that C v is 
numerically equal to k v . Similarly C v can be shown to be 
numerically equal to 2 k x . 

At low values of Z, C v approaches the value 2, which 
implies that k x =l and k v = 2. That these values of k rep- 
resent a critical combination of stresses for an infinitely 
long flat plate was shown in reference 12. At large values 
of Z, the curve approaches the curve given in figure 1 for 
buckling under lateral pressure alone and, like that curve, 
has branches representing buckling into two circumferential 
waves. 

In figure 5 the computed values of the pressure coefficient 
C p at which the cylinder would buckle if only the axial pres- 
sure were acting and if only lateral pressure were acting are 
compared with the results when both are acting because of 
hydrostatic pressure. At large values of Z the circumferen- 
tial stress at which buckling occurs under hydrostatic pres- 
sure is substantially the same as it would be if no axial stress 
were present, as in the case of lateral pressure. The reason 
that the circumferential stress appears as the main factor in 
buckling at high values of Z presumably is that at these 
values of Z the axial stress required to produce buckling is 
many times the circumferential stress required, whereas 
under hydrostatic pressure the axial stress actually present is 
only one-half the circumferential stress, 


prL z r _ prL 


r 


G 
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Figure 5. —Comparison of solution for buckling of simply supported cylinders under hydrostatic pressure with solutions for buckling under axial pressure alone and lateral 

pressure alone. 
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In figure 6 the curve of figure 4 is compared with curves 
representing Sturm's theoretical results (reference 5) and 
with a curve based on the following formula developed at 
the United States Experimental Model Basin (reference 13, 
equation (9)): 

r - ^ (ST y 

This formula is an approximation based on theoretical results 
obtained by Von Mises (reference 3, p. 479) which are iden- 
tical with the results in the present paper for buckling under 
hydrostatic pressure. Figure 6 shows that Sturm's theoret- 
ical results (reference 5) are in reasonable agreement with 
those of the present paper and that the formula from the 
United States Experimental Model Basin practically coincides 
with the present results except at very low values of Z. 

Test results from references 5 and 13 are included in fig- 
ure 6. The test data are in good agreement with the theoretical 
results except at low values of the curvature parameter Z 
at which the theoretical results are appreciably above those 
obtained experimentally. A possible explanation of the 
discrepancy between the theoretical and experimental results 
at low curvature is suggested by the relative importance of 
axial and circumferential stress in causing buckling. The 
axial stress becomes important only at low values of the 
curvature parameter Z. It is known experimentally that 
buckling under axial stresses may occur far below the theo- 
retical value of the critical stress. At low values of Z 


cylinders under hydrostatic pressure may therefore be 
expected to buckle well below the theoretical critical load 
just as cylinders do under axial compression. 

Torsion. — The problem of the determination of the buck- 
ling stresses of cylinders in torsion was solved by Donnell 
(reference 4) who gave an approximate solution of the equa- 
tion of equilibrium. A somewhat more accurate solution of 
this equation is given in reference 14. The essential results 
are shown in figure 7 taken from reference 14. At low 
values of Z the buckling-stress coefficient k s approaches the 
value 5.34 appropriate to infinitely long flat plates loaded in 
shear (reference 15). At higher values of Z the curve 
approaches a straight line given by 

k s =0.85Z 3/i 

At very high values of the curvature parameter the curve 
splits up into a number of other curves, depending on the 
value of r/t. The curves for various r/t values at high 
values of Z represent buckling into two circumferential 
waves. As mentioned before, Donnell's equation is not 

reliable for the case n = 2 ( a case which occurs for cylinders in 

torsion when >10 -0* A solution for this case given 

by Schwerin and discussed in reference 4 results in critical 
stresses about 20 percent below those of the present paper. 
Because Schwerin’s solution does not satisfy the condition 
w= 0 at the end of the cylinder, however, it is probable that 
the error in the present solution for n = 2 is less than 20 
percent. 



Fioure 6.— Comparison of present solution for buckling of simply supported cylinders under hydrostatic pressure with other theoretical solutions and test results. (Sturm’s results 

arc from reference 5 and Windenburg and Trilling’s results are from reference 13.) 
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In experimental investigations of cylinders in torsion the 
maximum rather than the critical loads have usually been 
reported. Because these maximum loads usually exceed the 
critical loads by only a small margin, it is common practice to 
check theoretical buckling stresses by comparison with the 
average stresses at maximum load. Such a comparison is 
provided in figure 8 which incorporates test data from refer- 
ences 4, 10, 16, and 17 For this figure the test results 
average about 15 percent below those given by theory. 

DISCUSSION 

Parameters appearing in buckling curves. The fact that 
the buckling of a cylinder under axial compression, lateral 
pressure, hydrostatic pressure, or torsion involves sub- 
stantially the same parameters is not a mere coincidence but 
is a direct consequence of the differential equation. The 
differential equation implies that when the requirement of 
an integral number of circumferential waves is removed the 
six variables L, r, t, E, n, and the load may be combined into 
two nondimensional parameters, one (k x , k v , k s , or C v ) 


describing the stress condition, and the other (Z) essentially 
determined by the geometry. (See appendix C.) It is also 
shown in appendix C that the buckling of a curved rec- 
tangular plate of any given length-width ratio may be 
represented in terms of these parameters. The critical 
stress of a cylinder or a curved plate of given length-width 
ratio may therefore be given by a single curve relating the 
two parameters provided that the number of circumferential 
waves may he regarded as continuously variable. This 
restriction becomes important at very large values of Z, for 
which the curves may split into a number of curves for 
cylinders of different values of r/t buckling into two circum- 
ferential waves. 

Except for hydrostatic pressure, each type of loading con- 
sidered results in a single uniform stress in the cylinder, and 
the nondimensional parameter k describing this stress is de- 
fined as follows in analogy to the parameter used in de- 
scribing the buckling of a flat plate: 

* (or r) 

tt 2 D 

LH 
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As the radius of the cylinder increases toward infinity (the 
other dimensions remaining constant), the cylinder approaches 
an infinitely long flat plate of the same thickness as the 
cylinder, having a width b equal to the length L of the 
cylinder. Accordingly, as the radius approaches infinity, 
the critical-stress coefficient k for the cylinder approaches the 
value of the corresponding stress coefficient for an infinitely 
long flat plate under the appropriate loading condition. 

The other nondimensional parameter Z is defined by the 
equation 

If the small correction due to Poisson’s ratio is neglected, a 
direct physical significance can be assigned to Z when its 
magnitude is small. The maximum distance from a slightly 
curved arc of length L and radius r to its chord can be shown 
to be given by the expression Z 2 /8r, which is called the" bulge” 
by some writers (see references 8 and 9). Accordingly, in 
the case of a curved strip of length L in the circumferential 


direction, L 2 /8rt is the bulge divided by the thickness and is 
thus a nondimensional measure of the deviation from flatness 
of the strip. As applied to a short cylinder, L 2 /8rt is the 
deviation from flatness of a square panel of the cylinder, 
each side of which is equal to the length of the cylinder. For 
cylinders having a length greater than a few tenths of the 
diameter, the parameter Z loses this simple physical signifi- 
cance and is perhaps best regarded as a nondimensional meas- 
ure of the length of the cylinder. Some indication of the 
variety of cylinder shapes corresponding to a fixed value of 
Z is given in figure 9. 

Boundary conditions. — When problems in the stability of 
cylindrical shells are solved by the use of DonneH’s equation, 
boundary conditions on u and v cannot be imposed directly 
because only w appears in the equation. The method of 
solution, however, may in some cases imply boundary condi- 
tions on u or v. In appendix D it is shown that for simply 
supported cylinders the method used in the present paper (a 
solution using one or more terms of a Fourier series satisfying 
the boundary conditions on w term by term) implies that at 
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— Constant L 


Figure 9.- Representative cylinders corresponding to the same value of 7 ( 7 . about 150). 


both ends of the cylinder the circumferential displacement i 
is zero, but that the cylinder edges are free to warp in the 
axial direction 0). For a simply supported rectangular 
curved panel, the present method implies (with regard to 
displacements within the panel median surface) zero displace- 
ment along the four edges of the panel and free warping nor- 
mal to the edges. These edge conditions on u and v are 
appropriate to cylinders or panels bounded by light bulk- 
heads or deep stiffeners which are stiff in their own planes 
but may be readily warped out of their planes. 

Relatively few calculations of the stability of a cylinder 
take into account the boundary conditions on u and v. A 
calculation for the case of torsion, however, was recently 
made by Leggett (reference 18). The results of this calcu- 
lation, computed for u=v=0 at the edges of the cylinder, are 
given only for Z<50. Throughout the range for which they 
are given, however, they agree very closely with the results 
found by the method employed in the present paper, which 
implies that at the edge of the cylinder z;=0 and 0. 
Restraining the ends of the cylinder from warping in the 
axial direction may therefore be assumed to have a negligible 
effect upon the buckling stress. This assumption receives 
added support from the form of the equation of equilibrium 
(appendix A) for the case of constant pressure 


l)\7 A w+t 



b 2 w , 9 b 2 w 
dx 2 “ T bxdy 


+ o* 


d 2 w 
by 2 + 


d 2 F 1 \ 
dx 2 r ) 


= 0 


In this equation, a x , (r y , and r are the stresses present just 
d 2 F . 

before buckling and >^ 2 is the circumferential stress produced 

by the buckling itself. The equation indicates that the only 
difference between the buckling behavior of a cylindrical 
sheet and that of a flat plate (found by omitting the last 
term in the foregoing equation) is due to the effect of the 
circumferential stresses caused by the buckling deformations. 


Because the restraint against warping in the axial direction 
requires the application of axial rather than circumferential 
stresses, this restraint might be expected to have only small 
effects on buckling stresses. Circumferential stresses would 
have to be applied to the straight sides of a curved strip to 
prevent warping normal to these edges during buckling. 
Because the circumferential stress due to buckling appears 
explicitly in the equation of equilibrium, the imposition of 
the restraint v = 0 to the straight sides of a panel should have 
an appreciable effect on the buckling stress (except when the 
straight sides of the panel are short compared with the curved 
sides). 

Theoretical results on the buckling of curved strips 
infinitely long in the axial direction are available to test the 
foregoing conclusion. In figure 10 the critical axial compres- 
sive stress for an infinitely long curved strip with u and v 
both zero along the edges (reference 8) is compared with the 
critical axial compressive stress when u is zero along the 
edges, and the edges are free to warp in the circumferential 
direction. (Sec appendix B for solution.) The critical axial 
stress is appreciably increased by the constraint v=0 in a 
certain range of small curvature. In figure 1 1 the critical 
shear stresses are compared under the same sets of edge 
conditions (references 6 and 7). The critical shear stress is 
conspicuously increased by the constraint p=0 except near 
the limiting case of flat plates. 

It appears from the foregoing discussion that the effect on 
the buckling stresses of preventing free warping normal to 
the curved edges of a cylinder or panel is very small but 
that the effect on the buckling stresses of a similar restraint 
on the straight edges of a panel may be quite important. 



Figure 10.— Comparison of the present, solution for the buckling under axial compression 
of a curved strip inflnitively long in the axial direction, with solution found by Leggett, 
(reference 8) . 
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Figure 11. — Comparison of Leggett’s solutions with present solutions for critical-shear-stress 
coefficients of a long curved strip. (Fig. 2 of reference 20.) 


Simplicity of results. — The theoretical results based on 
Donnell’s equation for the critical stresses of cylinders under 
a given loading condition appear particularly simple when 
presented as a logarithmic plot of buckling coefficient k 
against the curvature parameter Z. As r approaches 
infinity, and therefore as Z approaches zero, k approaches 
the value appropriate to a flat plate. At large values of Z 
the curve approached a straight line in each of the cases 
investigated. These straight lines had slopes 0.5, 0.75, and ] 
and are given approximately by the following equations 
which have already been given in the present paper and are 


reassembled here and provided with upper 
for easy reference: 

and lower limits 

k v = 1.04Z I/2 

(l00<z< j(t 


k r s =0.85Z 3/4 

(50 <Z< 10 (^ 


k x = 0.702Z 

(3 < z <s(l 

-Yd-,-)) 


These equations can also be written (when /z is taken to be 0.316) 


.,=0.926 f ( 5 )“- 

.0.926 25 (f )"(0 

0 

0 

*■*$ 

A 

f)< 5 0 

„ ... Et/tr\"‘ 
'=0.747 - (p) 

0 
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A 

9<“9 

Et 

<r x = 0.608 — 
r 


(a | < (#)< e 0 

II. MODIFIED EQUATION 

curvature. 

Equation (4) shows 


THEORY 

DERIVATION OF MODIFIED EQUATION 

The equation of equilibrium for a flat plate may be written 

( 2 ) 


nv 74 i ^ ( &W d 2 w , d 2 w\ 

DV w+t [a x dr ,+2r dxdy +« v dy ^P=° 


where j) is lateral pressure. (This equation is equivalent to 
equation (197) of reference 3.) 

For a cylindrically curved plate having a radius of curva- 
ture r, the following pair of simultaneous equations of equi- 
librium may be written (as a generalization of equations (11) 
and (10) of reference 7): 

nv74 . / d 2 w , d 2 w d 2 w\ 

DV w+t ^<r* bxi -r2r bxby + <ru by2 )+ 


t/VF, \ n 

P ~r\W* +,T '') =0 


V*F+-~ 2 = 0 
r ox 1 


(3) 


(4) 


where F is Airy’s stress function for the median-surface 
stresses produced by the buckle deformation (reference 19). 
Equation (3) differs from equation (2) only in the addition 
t /d 2 F \ 

of the term— — (-^2 d-o-j/b which expresses the effect of the 

823151 — 49 3 


cylindrical shells experience stretching of the median surface 
when originally straight lines in the surface are bent slightly. 
Elimination of F between equations (3) and (4) by suitable 
differentiations and additions gives the following single 
equation in w for the equilibrium of cylindrical shells: 


, Et d 4 w . . / 

Dv 8 w+~2 ^ +* v 


d 2 w _i_ 9 ^ 2w 

dx 2 7 dxdy 


+ 


d 2 W <Ty\ . 

r) + V P = () 


(5) 


Equation (5), which was first derived by Donnell (reference 4), 
was treated in part I. 

An alternative method for obtaining a single equation in w 
for the equilibrium of a cylindrical shell is to solve equa- 
tion (4) for F and substitute the result into equation (3). This 
procedure can readily be carried out in the following manner. 
Differentiation of equation (4) twice with respect to x gives 


4 d 2 F , E d 4 w 
dx 2 ^ r dx 4 


= 0 


( 6 ) 


d 2 F . 


The symbolic solution of equation (6) for is 

: _ A ' V 4 <>*w 


d 2 F 

dx 2 ' 


dx* 
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Substitution of this result into equation (3) gives 


DVho+ l 
r 



d 2 w 

dx 2 


+ 2r 


d 2 w 
dxdy ' 


d 2 w 
a “ dy 2 



(7) 


Equation (7) is simply equation (5) modified by multiplica- 
tion by the operator V -4 . In the present paper, equations (5) 
and (7) are referred to as Donnell’s equation and the modified 
equation, respectively. 


ADVANTAGES OF MODIFIED EQUATION 

One of the quickest and most convenient methods for ob- 
taining solutions of flat-plate buckling problems to any 
desired degree of approximation uses a Fourier series type of 
expansion for the deflection surface w. Both Donnell’s 
equation and the modified equation can be solved by this 
method in the case of buckling problems involving curved 
plates having simply supported edges. 

As mentioned in the “Introduction,” however, Donnell’s 
equation is not well adapted to solution by Fourier series of 
problems involving the stability of shells with clamped edges. 
The cause of the trouble appears to be that the calculation 
of some of the high-order derivatives found in Donnell’s 
equation sometimes leads to divergent trigonometric series 
when the edges are clamped. The modified equation, how- 
ever, is applicable to clamped-edge problems as well as to 
problems involving simply supported edges because lower- 
order derivatives are involved. 

Besides its advantages in the solution of problems involving 
shells with clamped edges, equation (7) has the additional 
advantage that each term has a definite physical significance: 
The first term gives the restoring force per unit area of the 
deflected surface due to bending and twisting stiffnesses; 
the second term gives the restoring force per unit area due 
to stretching stiffness; and the remaining terms give the 
deflecting or restoring forces per unit area due to applied loads. 
Because of these advantages, the modified equation was 
adopted for general use in references 11, 14, and 20 to 23. 

Both Donnell's equation and the modified equation result 
in the same critical stresses for simply supported cylindrical 
shells, and the two methods require essentially equivalent 
mathematical processes. (See appendix E.) The charac- 
teristics of solutions by means of Donnell's equation in the 
case of simply supported shells — namely, the theoretical 
cylinder parameters, the simplicity of calculations and re- 
sults, and the implied boundary conditions on u and v — are 
characteristics, also, of solutions by means of the modified 
equation. The same characteristics, except for a change in 
the implied boundary conditions on u and v , also apply to 
solutions of clamped-edge shell problems by means of the 
modified equation. This change is discussed in the section 
entitled “Boundary Conditions.” 


is that the Ritz method begins with an energy expression, 
whereas the Galerkin method begins with an equation of 
equilibrium. The Galerkin method is readily adaptable to 
the solution of equation (7) and is now described briefly. 

Let the equation of equilibrium be written 

<2(w) = 0 (8) 

where Q is some operator in x and y which for the purposes 
of this paper is taken to be linear. According to the Galerkin 
method, the equation may be solved by expanding the un- 
known function w in terms of a suitable set of functions 
Ji(x)gj(y), each of which satisfies the boundary conditions 
but not in general the equation of equilibrium: 

w=J2 X aufi(x)9j(y) (9) 

i j 

Substitution of this expression for w into equation (8) 
gives the following equation: 

0 ( 10 ) 

i j 

Because the functions Ji(x)gj(y) were chosen to satisfy 

the boundary conditions rather than the equation of equi- 
librium, equation (10) cannot, in general, be satisfied identi- 
cally by any choice of the coefficients a tj . These coefficients 
can be chosen, however, to assure the vanishing of certain 
weighted averages of the left-hand side of equation (10). 
The weighting functions used in the Galerkin method are 
the original expansion functions, so that the following 

simultaneous equations for determining the coefficients 
are obtained: 

2SS™#ii=0 (m= 1, 2, 3, . . .;n= 1,2,3,. . .) (11) 

i j 

where 

B mntj = ffMx)g n (y) Q[ft(x) gj {y)] dx dy (12) 

The simultaneous set of linear algebraic equations in the 
unknown coefficients a tj (equation(l 1 )), obtained by using 
the original expansion functions as weighting functions, is 
ordinarily the same set which would be found by the Ritz 
method, if the same series expansion for w were used. A 
solution of any desired degree of accuracy may therefore be 
obtained by the Galerkin method. 

In applying the Galerkin method to equation (7) by use 
of Fourier series expansion for w , expressions of the type 

_ . . IttX . jiry 

V - 4 XI X sm smY 

i j a 0 

must be evaluated. The operator V -4 , the inverse of V 4 , 
simply introduces into the denominator of each term of the 
series the expression that comes into the numerator if V 4 is 
applied. Thus, 


SOLUTION OF MODIFIED EQUATION BY GALERKIN METHOD 

An approximate method of solving vibration and buckling 
problems closely paralleling that of Ritz was introduced in 
1915 by Galerkin. (See, for example, references 24 and 25.) 
The main distinction between the Ritz and Galerkin methods 


v 4 X X a u sin — sin -T- = 

it a b 


XX 




. ITX • 77T?/ 

sin sm — 
a b 


(13) 
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This result may readily be verified by applying the operator 
V 4 to each side of equation (13). 

In writing equation (13) the quantity V -4 /, as defined by 
the equation 

v 4 v- 4 /=/ 


was tacitly assumed to be unique. The quantity actually 
is not unique; any number of terms which vanish when 
operated upon by V 4 could be added to the right-hand side 
of equation (13). The omission of such terms makes the 
present analysis parallel to the analysis using Donnell’s 
equation (see part I) and implies certain boundary conditions 
on u and v, which are discussed in a subsequent section 
entitled “Boundary Conditions.” 

DEFLECTION FUNCTIONS 

Simply supported edges. — For simply supported cylindri- 
cal shells, the following series expansions for w may be used 
to represent the buckle deformation to any desired degree of 
accuracy (in these functions, x is the coordinate in the axial 
direction and y , the coordinate in the circumferential 
direction): 

(1) Rectangular curved plate (axial dimension a and 
circumferential dimension b) 


w= T! a 

TTi—l n — 


. mirx . nivv 
mn sin — - sin — 


(14) 


(2) Curved strip long in the axial direction (circumferen- 
tial width b and axial wave length 2\) 

(a) Direct stresses only 


. irx ^ . miry 

w = sin a m sm — - 

A m = 1 


T- < i5 > 

(b) Shear stress with or without addition of direct stress 


TTX 


rrnry 


TTX 


miry 


. n ju . non if . n Us 7 . ni/liu 

w= sm 2-j dm sm - ; * +cos -r 2 j o m sm - , (16) 

A 771=1 o A 7 7i=l 0 


(3) Complete cylinder (length L and circumferential wave 
length 2X) 

(a) Direct stresses only 

. mi . mirx 

W=sm 2D d m sm -j- (i 7 ) 

A 771=1 -L v ' 

(b) Shear stress with or without addition of direct stress 


w = sm 


iry 


dm sm 


mirx 


Try \ 7 
COS 2_J On 
A m = l 


. mirx . . 

sm— j— (18) 


Clamped edges. -Probably the simplest method of treating 
cylindrical shells with clamped edges is to employ the 
expansions in equations (14) to (18) modified by substituting 
functions of the type 


, \ . mirx . irx 1 f (m-l)irx (m+ lWafl 

^W=sm — sm — =2 [cos cos - v — ' J 

(19) 

wherever functions of the type sin * 1 ^ appear, with a 

similar substitution for functions of y (all terms involving 
summation subscripts m and n are thus changed; terms 


involving X, such as sin remain unchanged). The func- 

tions <p m (x) form a complete set so that finite expansions for 
w of the type suggested for shells with clamped edges as well 
as those for shells with simply supported edges may be used 
to represent the buckle deformation to any desired degree of 
accuracy. 

BOUNDARY CONDITIONS 


Simply supported edges. — Appendix D shows that, if the 
buckling stress of a simply supported shell is found by means 
of the expansions for w given in the preceding section en- 
titled “Deflection Function,” the boundary or edge condi- 
tions implied for the median-surface displacements u and v 
are zero displacement along each of the edges of a cylinder 
or curved panel and free displacement normal to each edge. 
(Although the proof given used equation (5), the proof 
could equally well have been based on equation (7).) 

The boundary conditions for simple support may thus be 
written, at a curved edge (x — Constant), 


d 2 w d 2 F 


( 20 ) 


and, at a straight edge (y — Constant), 


d 2 w d 2 F 
w =dy* =U =dx> =0 


( 21 ) 


Clamped edges. — By a method similar to that in appendix D 
solutions using the functions suggested in the preceding 
section for the treatment of clamped edges can be shown to 
correspond to the boundary conditions zero displacement 
normal to an edge and free displacement along an edge. 

The boundary conditions for clamped edges therefore be- 
come, at a curved edge (x= Constant), 


dw d 2 F 

w= ^ =u= ^ 2=0 


( 22 ) 


and, at a straight edge (y— Constant), 

dw d 2 F / OQ \ 

Wz= ^ =V= d^ = 0 (23) 

Discussion. — As mentioned in part I, the boundary con- 
ditions implied for u and v in the case of simply supported 
edges are appropriate for cylinders or panels bounded by light 
bulkheads or deep stiffeners, which are stiff in their own 
planes but may be readily warped out of their planes. 

The boundary conditions on u and v appropriate for a 
clamped edge would seem to be zero displacement normal 
to the edge and zero, rather than free, displacement along 
the edge. Comparison of critical stresses for shells with 
clamped edges found by the method in the present paper 
with critical stresses found by the method in references 7 
and 8, giving boundary conditions u=v= 0, however, indi- 
cates that the imposition of the added requirement of zero 
displacement along the edge ordinarily has very little effect 
on the critical stresses. 

A less satisfactory method for solving problems concerning 
shells with clamped edges involves the use of functions of the 
type 


1 . mirx 

sm 
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sm 
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instead of those described by equation (19). In this method, 
the functions used are those for simple support taken in 
such combinations that the edge slope is zero. Use of such 
functions leads to the same boundary conditions on u and 
v as were described for simply supported edges; at the edge 
y= Constant, for example, the boundary conditions become 


dw d-F 

=u= ^ 2 = 0 

d?/ ox z 


(24) 


The use of these functions to represent shells with (damped 
edges is not recommended, however, for the following reasons: 
The associated boundary conditions seem to be artificial 
and unlikely to be reproduced even approximately in actual 
construction; the method leads in some cases to solutions that 
differ considerably from the solution for ideal elamped-edge 
conditions in which u=v=0; and the solutions obtained 
generally converge rather poorly. 

APPLICATIONS AND DISCUSSION 

Among the more difficult shell-stability problems to treat 
theoretically are those which involve shear stresses. In fact, 
until 1934, when Donnell's paper on critical shear stress of 
a cylinder in torsion was published (reference 4), such prob- 
lems were generally regarded as impracticable to solve. In 
order to illustrate the type of solution to be found by the 
method of analysis just outlined and the effect of boundary 
conditions on critical stresses, the results obtained for a 
number of shell-stability problems involving shear stresses 
are reproduced and discussed briefly here. The problems 
treated are summarized in table I. 

Critical shear stress of long curved strip. — The critical 
shear stress for a long plate with transverse curvature is given 
by the equation 


where k s is a dimensionless coefficient, the value of which 
depends upon the dimensions of the strip, Poisson's ratio 
for the material, and the type of edge support. In figure 12 
(fig. 1 of reference 20) the shear-stress coefficient k s is given 
for plates with simply supported edges and with clamped 
edges. This solution for simply supported edges coincides 
with that given by Kromm (reference 6). 

As indicated in the previous section entitled “ Boundary 
Conditions," the solution corresponding to the boundary 
conditions of equation (24) (dashed curve of fig. 12) is poorly 
convergent and deviates appreciably from the results for 
completely fixed edges. Figure 12 shows this poor conver- 
gence in the limiting case of a flat plate, for which the critical 
stress is independent of boundary conditions on u and v. 
Even a tenth-order determinant led to a result that is 
7 percent above the true solution; whereas the result using a 
fourth-order determinant obtained with the deflection 
functions recommended for clamped edges is only 1 percent 
above. 

In figure 11 (fig. 2 of reference 20) the solutions given in 
figure 12 are compared with the results given by Leggett 
(reference 7) for simply supported and clamped edges with 
u =v = 0 at each edge. Throughout the range for which 


TABLE I— INDEX OF PROBLEMS TREATED 


Problem 


y 


v_y 

'Qi 

jpr 

j i 

n 

Q 

w 


Figure 

Reference 

12 

20 

11 

20 


7 

(LeggetO 

7 

14 

13 

21 

14 

21 

Not 

shown 

22 

15 (a) 

23 

15 (b) 

23 


(Simply supported (u- 0, v^O). 
^Clamped (u^O, v=0). 
(.Clamped (m = 0, v^O). 


Simply supported (?/ = 0, v^O). 
Clamped (u^O. v=0). 
Clamped (u=0, v^O). 

(Simply supported (?/ = t>=0). 
\Clamped (u=v= 0). 


Edge condition 


(Simply supported. 
\ Clamped. 


Simply supported. 


Simply supported. 


(Simply supported. 
\ Clamped. 


Simply supported. 


Clamped. 


they are given, Leggett's results for clamped edges differ 
only slightly from those of the present paper. On the other 
hand, the previously mentioned discrepancy between the 
results for completely fixed edges (?/ = y=0) and those for 
the boundary conditions of equation (24) (dashed curve) 
may be inferred from this figure to be considerable for large 
values of Z. A minimum measure of this discrepancy is the 
distance between the clamped-edge curves for and for 
u = 0 in figure 11, since Leggett's curve must always lie 
above the curve for v=0. 

The reason for the marked increase in buckling stress of 
simply supported curved strips when the edges are restrained 
against circumferential displacement during buckling is 
discussed in part I. 

Critical shear stress of cylinder in torsion. -The critical 
shear stress of a cylinder subjected to torsion is given by the 
equation 



In figure 7 (fig. 1 of reference 14) the values of k s are given 
for cylinders with simply supported edges (boundary con- 
ditions of equation (20)) and cylinders with clamped edges 
(boundary conditions of equation (22)). At high values 
of Z, the values of k s for thick cylinders are given by special 

curves for various values of y -y/l — m 2 , as discussed in part I. 

At values of Z greater than about 100 only a small 
increase in buckling stress is caused by clamping the edges. 
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Fic.uke 1 2 . —Critical-shear-stress coefficients for a long curved strip. (Fig. 1 of reference 20.) 


The results indicated in figure 7 are in very close agreement 
with Donnell's results for the same problem, except in the 
range 5 500 where the somewhat lower curves of the 

present paper represent a more accurate solution. 

Part 1 shows that boundary conditions imposed upon u 
and v at the curved edges of a panel or cylinder have an 
almost insignificant effect on the buckling stresses, whereas 
conditions imposed on v at the straight edges may be quite 
important. Comparison of figure 12, in which boundary 


conditions on straight edges are considered, with figure 7, in 
which conditions on curved edges are considered, indicates 
that a similar situation exists with respect to restraint 
against edge rotation. 

Critical shear stress of curved panel. — The values of k s 
giving the critical shear stresses of simply supported curved 
rectangular panels are given in figures 13 and 14 (figs. 1 and 
2, respectively, of reference 21). The corresponding bound- 
ary conditions on u and v are zero displacement parallel to 



Figure 13.— Critical-shear-stress coefficients for simply supported curved panels having cir 
cumferential dimension greater than axial dimension. (Fig. 1 of reference 21.) 



Figure 14.— Critical-shear-stress coefficients of simply supported curved panels having axial 
dimensions greater than circumferential dimension. (Dashed curve estimated.) (Fig. 2 
of reference 21.) 
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the edges and free warping normal to the edges. Figure 13 
indicates that, as the curvature parameter Z increases, the 
critical shear stresses of panels having a circumferential 
dimension greater than the axial dimension approach those 
for a cylinder. Figure 14 indicates that, as the curvature 
parameter Z increases, the critical shear stresses for panels 
having an axial dimension greater than the circumferential 
dimension deviate more and more from the critical shear 
stress for an infinitely long curved plate. Reference 21 
shows that the reason for this deviation in figure 14 is that 
at high curvatures the buckling stresses of these panels, as 
well as those of figure 13, approach those of the cylinder 
obtained by extending the circumferential dimensions ol the 
panels. 

The effects of boundary conditions in the limiting cases of 
infinitely long curved strips (fig. 12) and of complete cylin- 
ders (fig. 7) suggest that the curves of figure 13 are sub- 
stantially independent of edge restraint at large values ol Z 
but that the curves of figure 14 would be considerably 
affected by a change in edge restraint. 

Long curved strips under combined shear and direct axial 
stress. Reference 22 shows that the theoretical interaction 
curve for a long curved strip under combined shear stress 
and direct axial stress is approximately parabolic when the 
edges are either simply supported or clamped, regardless of 
the value of Z. This parabola is given by the formula 

Rs 2 +R x = l 

where R s and R x are the shear-stress and compressive-stress 
ratios, respectively. 

At high values of Z curved strips, like cylinders, buckle at 
compressive stresses considerably below the theoretical crit- 
ical stresses. In order to take this condition into account, 
certain modifications in the theoretical results are proposed 
in reference 22 for use in design. 

Cylinders under combined shear and direct axial stress. 
The theoretically determined combinations of shear stress 
and direct axial stress which cause a cylinder with simply 
supported and clamped edges to buckle are shown in figure 
15 (fig. 1 of reference 23). Considerable variation in the 
shape of the interaction curves occurs for low values of Z. 
For high values of Z the interaction curves for either simply 
supported or clamped edges are similar to the curve for 
Z=30. 

Because cylinders actually buckle at a small fraction of 
their theoretical critical compressive stress, the theoretical 
interaction curves of figure 15 cannot be expected to be in 
satisfactory agreement with experiment when a very ap- 
preciable amount of compression is present. For semi- 
empirical curves and a check of available test data, see 
reference 23. 

CONCLUDING REMARKS 

The use of Donnell’s equation to find the buckling stresses 
of simply supported cylindrical shells leads to simpler results 
and involves less labor than the use ol equations in which 
second-order terms are retained. The buckling stresses 
found by use of Donnell’s equation are in reasonable agree- 
ment with results based on other theoretical calculations. 



Figure is.— C ritical combinations of shear-stress and direct-axial-stress coefficients for 
cylinders. (Fig. l of reference 23.) 


Except for the case of axial loading, they are also in reason- 
able agreement with test results. Boundary conditions having 
to do with axial and circumferential displacements cannot 
be handled directly by use* of Donnell’s equation. This 
disadvantage is not considered serious, however, because the 
boundary conditions on axial and circumferential displace- 
ment, which are implied by the simple solutions given, 
correspond approximately to those that are most likely to 
occur in practical construction and because in many cases 
the buckling stress is not very sensitive to these boundary 
conditions. The restriction to simply supported edges in 
Donnell’s equation can be removed by the introduction ol 
a new equation which is equivalent to Donnell’s equation 
but is better adapted to solution by Fourier series. This 
modified equation can be solved for the buckling stresses of 
curved sheet having either simply supported or clamped 
edges by established methods essentially equivalent to those 
in use for fiat sheet. This approach permits a simple and 
straightforward solution to be given for a number of prob- 
lems previously considered rather formidable. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 20 , 1947. 


APPENDIX A 

SIMPLIFIED EQUATIONS OF EQUILIBRIUM FOR CYLINDRICAL SHELLS 


The principal sets of simplified equations currently in use 
for the equilibrium of cylindrical shells are listed for con- 
venient reference. The various sets of equations are equiva- 
lent. The reference papers in which the equations are derived 
are also listed. The equations given arc generally not 
identical with those in the reference papers but are modified 
in certain respects to include all the loading conditions 
studied in the present paper or to put them in the notation of 
the present paper. 

The three following simultaneous equations in displace- 
ments u, v, and w (reference 3) are derived from the conditions 
of static equilibrium: 


c )Hl 1 —y.d 2 u . 1+m /x dw_ 

dx 2 ' 2 dy 2 ' 2 dxdy^~ r dx 

d 2 v 1— At d 2 v , 1+m d 2 u , 1 dw 
dy 2 2 dx 2 2 dxdy r dy 


. JSi /do , du w\ 

Dv ”+Ki-r)U + "^ + <) + 


/ d 2 w. n d 2 w d 2 w a- A, 

\ ' r dx 2 ' wT dxdy ' (Jy dy 2 r ) ^ 


9 = 0 


(Al) 
(A 2) 

(A3) 


Two simultaneous equations in deflection w and stress func- 
tion F (reference 6) are as follows: 


d 2 w 


r dx 2 


d 2 w 


(A 4) 


UV'w+1 a x d/ , + 2 t dj , dy + Gy dy , r df2 + 


V = 0 (A 5) 

A single equation in deflection w (Donnell’s equation, refer- 
ence 4) is 

d 2 w . d 2 w 


T-\_n . Et d'w , < / 


(Tx dx 2 ^ “ 7 dxdy 


d 2 w o \ 

7 by 2 


T /)+Vp=0 (AO) 


The relationships between u and w and between v and w art* 
(reference 4) 


1 d 3 w 
'dxdy 2 


(A7) 

v d 3 w 

d*w 

(AS) 

' dx 2 dy 
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APPENDIX B 


Donnell's equation for the equilibrium of cylindrical shells 
is used to investigate the stability of simply supported 
cylinders subject to lateral pressure, axial compression, and 
hydrostatic pressure, and of simply supported curved strips 
long in the axial direction subject to axial compression. 

CYLINDER UNDER LATERAL PRESSURE 

If bending of the cylinder wall is neglected, constant 
lateral pressure on a cylinder causes only circumferential 
stresses. Donnell's equation (equation (AO)) then reduces 
to 


THEORETICAL SOLUTIONS 

where 


where 


7i _ s , Et w . ,_ 4 dhv _ 

DVw +r 2 dif 


„ -Vr 
**- t 


(Bl) 


and p is the pressure applied. (By virtue of the preceding 

equation the terms involving p and !/ appearing in equation 

(A6) cancel in the present case.) Division of equation (Bl) 
by D results, with proper substitutions, in the following 
equation: 


8 12 Z 2 d A w , 7 r 2 4 d 2 w 

V w+- l t- d -4 +K - D v d - 2 - 0 


(B2) 


The boundary conditions corresponding to simply supported 
edges (no deflection and no moment along the edges) arc 

w(0,y)=w(L,y)=0 

!?( O.V)=5?(*,V>=0 

A solution of equation (B2) satisfying t lie boundary condi- 
tions for simple support is 


. 7T?/ . m TTX 

w=w 0 sm 1 sm r 

A 1 j 


(B3) 


where X is the half wave length in the circumferential direc- 
tion. Combining equation (B3) and equation (B2) yields 
the following equation: 


(m 2 + /3 2 ) 4 + 1 2Z l m, -k u F(m 2 +py=() 

7T 

The solution of equation (B4) for k y is 

7 _ (m 2 + d 2 ) 2 | 12Z 2 ?n 4 

P 2 7 r 4 /3 2 (m 2 +/3 2 ) 2 

IS 


(B4) 


(B5) 


0= 


L 


The critical value for k y is found by minimizing the right- 
hand side of equation (B5) with respect to m and 0. If the 
numerator and denominator of the last term in equation 
(B5) are divided by m 4 , it becomes evident that under the 
restriction of integral values of m , k v will be a minimum 
when m— 1. Equation (B5) therefore becomes 


h s i +py i 


12Z 2 

,2 (l+d 2 ) 2 


(BO) 


The results found by minimizing this expression for k y with 
respect to (considered continuously variable) are shown in 
figure 1 by the curve independent of r/t. 

At low values of Z, buckling is characterized by a large 
number of circumferential waves. As Z increases, the num- 
ber of circumferential waves decreases until it finally becomes 

two (*-5> corresponding to buckling into an elliptical 

cross section. The curves for buckling into two circumferen- 
tial waves are shown in figure 1 as the curves for various 

values of y -y/ 1 — /x 2 . The equations for these curves arc 

found by substituting in equation (B5) the last of the follow- 
ing expressions for 

L = 2L = 2 I Z 

CYLINDER IN AXIAL COMPRESSION 

When only axial stress is present, equation (AG) becomes 


, Et d 4 w . . . d 2 w 

Dv + dx 2 = 0 


(B7) 


Division by D results, with proper substitutions, in the 
follow r ing equation: 


„ \2Z' 2 b i w. J ir 2 _.d 2 w 

v»+ s 4i,„v d ,.-o 


<B8) 


Combination of the deflection equation (B3) with equation 
(B8) yields the following equation: 

(m 2 +f} 2 y+ l - Z E-k x m 2 (m 2 +l } 2 ) 2 = 0 (B9) 

7T 
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The solution of equation (B9) for k x is 

(m 2 +0 2 ) 2 12Z 2 m 2 

kx ~ m 2 r 4 (m 2 +0 2 ) 2 

The critical value of k x for a given value of Z may be found 
by minimizing k x with respect to the parameter 

(m 2 +0 2 ) 2 

m 2 

If no restrictions are placed on the value that this parameter 
can take, the minimum value of k x is found to be 


^ ( m 2 -}~0 2 ) 2 12Z 2 m 4 


(BI3) 


The critical value of C v is found by minimizing the right- 
hand side of equation (B13) with respect to m and 0, with 
due regard to the values which rn and 0 may assume. It 
can be shown that the minimum value of C v is found by 
taking m equal to 1, so that equation (B13) becomes 


r _(l + 0 2 ) 2 12Z 2 

1 • /I \ 

i+/3 2 tt 4 (1 + 0 2 ) 2 (i + 0 2 J 


(B14) 


*,=-■>4- Z=0.702Z (BlO) 

7 T“ 

which coincides with the results generally given for the 
buckling of long cylinders. 

For values of Z below 2.85, however, the straight-line 
formula (equation (BlO)) cannot be used, since it implies 
either imaginary values of the circumferential wave length X 
or the number of axial half waves m below unity. The 
critical stress coefficient k x for Z<2.85 is found by substitut- 
ing the limiting values 0=0 and m= 1 in equation (B9). 
The results are shown in figure 3. 

CYLINDER UNDER HYDROSTATIC PRESSURE 

Hydrostatic pressure applied to a closed cylinder produces 
the following axial and circumferential stresses: 




pr 

T 


The. equation of equilibrium (equation (A6)) when both 
circumferential and axial stress are present is (since V 4 p = 0) 


. Et d 4 w , ^ d 2 w . . d 2 w „ /T -. * 

DV <r x tV y c 2 “l“ 0 V^ 7 5y 2 ^ (BH) 


r 2 dx 

t 

By use of the definition 


r —VrL 2 

'Dir 2 


equation (Bll) can be written 
vw + IS U 




If the deflection equation (equation (B3)) is combined with 
equation (Bl2), the following expression results for C p : 


^Equation (B14) is equivalent to an equation derived by 
Von Mises (reference 3, p. 479). The results of minimizing 
V p with respect to 0 are shown in figure 4. (The curves 

given for various values of -- V 1 — M 2 have the same signif- 
icance as in the case of a cylinder buckling under lateral 
pressure alone.) 

LONG CURVED STRIP IN AXIAL COMPRESSION 

Because it merely describes equilibrium at a point, equa- 
tion (B7) applies to the buckling of a long curved strip as 
well as to cylinder buckling. In modifying this equation to 
obtain nondimensional coefficients as in equation (B8), 
however, it is convenient to define k x and Z in terms of the 
width of the strip b rather than in terms of the axial length 
L , which applied in the case of the cylinder. Accordingly, 
equations (B7) and (B8) for a cylinder in axial compression 
may be applied also to the buckling of a curved strip, long 
in the axial direction, subjected to axial compression, 
provided the curved width b is everywhere substituted 
for the length L. Substitution of the deflection 


. 7 tx . mry 

W=w 0 Sill y Sin — y 


into equation (B2) (modified by substitution of b for L) gives 


k x = 


(n 2 +0 2 ) 2 , 12 Z 2 0 2 


0 2 


where 


V+0 2 ) 2 


(BIS) 



Equation (B15) is very similar to equation (B9) and each 
equation yields the same critical value for k x at large values 
of Z. At small values of Z, the minimum value of k x is 
found by taking n= 1 in equation (B15) and minimizing 
with respect to 0 the result ing expression for k x . The results 
are given in figure 10 together with results found by Leggett 
(reference 8). 


APPENDIX C 


PARAMETERS 


It is shown that Donnell's equation implies that under 
certain limitations the buckling coefficient k , familiar from 
flat-plate theory, can be expressed in terms of the curvature 
parameter Z alone in the case of a complete cylinder or a 
curved rectangular panel of given length-width ratio. 

DonnelFs equation (A6) is (when j) is constant or zero) 


Even without solving this equation it is clear that w must 
be a function of the independent variables £ and rj, and also 
the parameters Z , k x , k s , and k v , and the derivatives of w 
will be functions of the same variables and parameters. 
Thus, if only one type of loading (represented by the buckling 
coefficient k) is present, equation (C2) may be written 


t ~ v_^o . Etd 4 w . / d 2 w . d 2 w . d 2 w\ „ , r 

DV W+ r 2 +N V* dx- + 2t dxdy + ( ° 0 


Let 


6 « 


y 

b =r] 


and 


Then 


v 2 _ a 2 

Vo -be + drj* 


V 0 2 =6 


\dx-^dy-J 


Multiplication of equation (Cl) by b 8 and substitution of 
the dimensionless coordinates £ and rj gives 


o , Etb 4 d 4 w . . / d 2 w 

DV C S W + )2 +b-tV ( / { rT x 

Division by D results in 


d 2 w 

d&y 


„ — . d 2 «A 


: 0 


8 Etb 4 d 4 w b 2 t 4 / d 2 ic 5 2 uA 

Vc w+ D? W + D V ° v x d£ 2+2r d£d„ + <r " dr' ) = 


JS? 


or, since D = , 0 " 

12(1— m ) 

x-7 8 I 1079^ I -n 4 // 1 0 / 1 / ^ 2 w\ , 

V 0 W+12Z- J^ + T-Vo (*« d?2 a,* )=< 


where 


(C2) 


7 _& 2 , , 

Z ~rt V1_M 


fi(£,y,Z,k) + \2Z 2 f 2 ^,rj,Z,k)-]-Tr 2 kJ^,r},Z,k) — 0 (C3) 

where and / 3 are definite, though unknown, functions. 
The variables £ and rj may now be eliminated by integration 
of both sides of this equation over the entire range of £ and 77. 
In the case of a curved panel of circumferential dimension a 
and axial dimension b the resulting equation is 

a 

* 2 kMl,v,Z,k)] = 0 (C4) 

The integrals of the functions f u f 2l and / 3 depend only 
upon Z, k, and the value of the ratio a/6. Accordingly, 
equation (C4) implies that a relationship of the following 
type exists: 

jU(k,zfy = 0 (C5) 

Equation (C5) indicates that for any given value of the panel 
aspect ratio a/6, the critical-stress coefficient k depends only 
upon Z. 

If a complete cylinder of length L rather than a panel of 
length 6 is under consideration, and t lie deflection w is 
periodic with wave length 2\ in the circumferential coor- 
dinate, the integration 

a 

rb 

1 ,h 

appearing in equation (C4) may be replaced by 

2\ 

J> 



where £ and 77 are now defined as xjL and yjL, respectively. 
The result then becomes 


Dt 2 

__ <r„tb 2 
v ~ Dt 2 


/. (k,z, 2 £)=o 

k = h(z, (C6) 
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The actual buckling stress is found by minimizing k with 
respect to 2 \/L. 

Theoretically, X must satisfy the equation 

irr=n\ (C7) 

where n is the number of circumferential waves and therefore 
an integer. When many circumferential waves are present, 
however, this restriction does not significantly affect the 
buckling stress, and the minimization of k with respect to 
2\ 

(considered continuously variable) leads to the result 
b 

k=MZ ) (C8) 

Equation (C8) indicates that provided the number of cir- 
cumferential waves is not too small the critical-stress coeffi- 
cient for a cylinder depends for practical purposes only upon 
the curvature parameter Z. 


When n is so small that its integral character must be taken 
into account, it appears from equations (C6) and (C7) that 
k depends upon both Z and r/X. Since, however, 



k for small values of n can alternatively be expressed in terms 

y . 

of Z and -- yl —jjl 2 , as in figures 1, 4, and 7. 
z 

By a similar analysis, it can be shown that when the buck- 
ling of a cylinder under hydrostatic pressure is represented 
by plotting the pressure coefficient C P against Z, a single 
curve is obtained except where the small number of circum- 
ferential waves requires splitting the curve into a series of 

curves for different values of — yl— y, 2 . 


APPENDIX D 


BOUNDARY CONDITIONS ON EDGE DISPLACEMENTS WITHIN THE MEDIAN SURFACE 


The solution of Donnell’s eighth-order partial differential 
equation for the stability of cylindrical shells is not unique 
under the imposition of the ordinary boundary conditions 
for simply supported or clamped edges. Two more bound- 
ary conditions at each edge, for example, one condition for 
u and one for v, are required to define completely the phys- 
ical problem and are therefore needed to make the solution 
unique. Because only w appears in the equation, boundary 
conditions on u and v cannot be imposed directly; they may, 
however, be implied by the method of solution. The pur- 
pose of this appendix is to show what boundary conditions 
on u and v are implied by the method of solution used in 
the present paper. In order to simplify the discussion, the 
analysis will first be made for the case w r hen only axial 
compression is present and will then be extended to other 
cases. 

When only axial stress is present, Donnell’s equation 
(equation (A6)) becomes 


DV*w+ 


Et d A w 
r l dz 4 


+ <t^V 4 


d 2 w 

dx 2 


= 0 


If the shell described by this equation is a curved panel with 
the origin of coordinates in one corner of the panel, a solu- 
tion satisfying the usual boundary conditions for simple 
support is 


. mwx . utv 
w=Wq sin - sin 

a b 


(Dl) 


where m and n are integers. This solution is also the 
solution to the problem of the buckling of an infinite two- 
dimensional array of panels identical to the one under con- 
sideration. (See fig. 16.) When such an array buckles, the 
displacements u, v, and w as v eil as the stresses, described 


by the stress function F, may he presumed to he periodic 
over the interval 2 a in the axial direction and 2b in the 
circumferential direction. 

Any function u(x, y) that is periodic with a wave length 
2 a in the ^-direction and with a wave length 2b in the 
^-direction may be expanded as follows (see, for example, 
reference 26) : 


U= ]C a„ 


. rmrx . utv , 
sm a. sm b + 


00 00 


S S b,„n sin 

VI = 1 72 = 0 


rmrx utv , 

cos — j- -- + 
a b 


00 00 



in = 0 7i = l 




ran 


7/2 = 0 71=0 


TrnrX . UtV . 
COS Sin “j~ 

a b 


m tx mry 

cos cos j ' 
a b 


(D2) 
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The relationship 
(equation (A7)) 


which must exist between u and w is 


rV*u= — ix 


dx 3 + 


d z u> 

dxdy 2 


( D11 ) 

(a,y) = 0 (D12) 


Substitution into this equation of the expressions for u and w 
from equations (D2) and (Dl), respectively, and use of the 
orthogonality of the functions in equation (D2) leads to the 
result 



/ TUt\ 
l \a ) 


ill 

= “'1 

[(T)‘+(t)J 


m 7r.r 

a 


sin 


n iry 

b 


Accordingly, the boundary conditions on u a: 

re 

u(x,0) = 0 

(D3) 

u(x,b)= 0 

(D4) 


(D5) 


(D6) 

Similarly by use of equation (AS) instead ol 
it can be shown that the boundary conditions 

' equation (A7) 
on v are 

v(0,y ) =0 

(D7) 

v(a,y)=0 

(DS) 


(DO) 


(DIO) 


The boundary conditions of equations (D5), (D6), (DO), 
and (DIO) may be combined to give four boundary condi- 
tions on the stresses induced by buckling. These boundary 
conditions, which are also derivable from equation (A4) by 
a method analogous to that just used to derive the conditions 
relating to u, are 


b 2 F 


(*, 0)=0 

b 2 F 

(X,b)= 0 


(D13) 


(D14) 


i b 2 F . c ) 2 F 
where A > and A 2 are, 
by- ox 2 


respectively, the median-surface 


axial and circumferential stresses caused by buckling. The 
eight boundary conditions given by equations (D3), (D4), 
(D7), (D8), and equations (Dll) to (D14), plus the eight 
boundary conditions on w for simple support of the four panel 
edges taken together uniquely determine the buckling stress. 

Although the preceding discussion of boundary conditions 
started with the assumption of axial stress only, the only use 
made of this assumption was in obtaining equation (Dl) as 
the solution for the buckling deformation. The same defor- 
mation, and hence the same arguments, apply when circum- 
ferential stress is present. When shear is present, a series of 
terms of the type in equation (Dl) must be used to repre- 
sent the deflection surface, and hence series of terms occur 
in the expressions for u, v, and F. Since the boundary condi- 
tions derived in the preceding analysis apply to each of the 
terms individually, by the principle of superposition they 
must also apply for the sum, so that equations (Dll) to 
(1)14) represent the boundary condition no matter what the 
applied stresses are. 

In summary it may be stated thatthe substitution of one 
or more terms of a double-sine-series expansion for w into 
Donnell’s equation and solution of the resulting equation for 
the buckling stress gives the solution corresponding to the 
following boundary conditions: 

(1) Each edge of the panel (or cylinder) is simply sup- 
ported; that is, the displacement normal to the surface of the 
panel and the applied moments are zero at the edges. 

(2) Motion parallel to each edge during buckling is 
prevented entirely. 

(3) Motion normal to each edge in the plane of the sheet 
occurs freelv. 


APPENDIX E 


COMPARISON OF RESULTS OBTAINED BY USING DONNELL’S EQUATION AND THE MODIFIED EQUATION IN THE 

STABILITY ANALYSIS OF SIMPLY SUPPORTED CURVED PANELS 


SOLUTION OF DONNELL’S EQUATION 

Donnell’s equation expressing the equilibrium of a curved 
panel under constant median-surface stresses can be written 
in general form as 

. Etb*w . , „ rot d 2 w _ 4 d 2 w_ . 

DVw + r 2 fai + vJV 5 x 2+2t<V + dy , 0 (El) 

where x is the axial coordinate and y the circumferential 
coordinate. Division of equation (El) by D and the intro- 
duction of the dimensionless stress coefficients k x , k y , and 
k s , and the curvature parameter Z results in the following 
equation: 


_ ,12 Z 2 d%. 7 7T 2 T74 d 2 ^, oZ 7T 2 4 p 2 W , 7 7T 2 _ 74 d 2 ^_ n 
VwH p- dx 4 +k x b *V dx 2 + 2 ^ 6 2 V dxby +k ' J b 2 ^ dy 2 ~° 

(E2) 


diere 
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k —a 

Kx—Vx 

k~r- bH 
ICS ~ T TT 2 D 

l _ bH 

kv-Vv t 2 D 


7 'i 2 

Z= H Vl_/1 


Equation (E2) can he represented by 

<2!(w)=0 

where is defined as the operator 


(E3) 


IlD- 0‘ ^ d 2 , r 2 y 2L 2 V 4 y 

V + b* dx i+kx b 2 ' dx 2± ^ s b 2 dxd>y " 6- D?/ 2 

The equation of equilibrium (equation (E3)) is solved by 
using the Galerkin method as described in the section 
entitled “Theory.” In applying this method the unknown 
deflection w is represented in terms of a set of functions (see 
equation (9)), each of which satisfies the boundary conditions 
but not in general the equation of equilibrium. A suitable 


set of functions of this type, which satisfies the boundary 
conditions for simple support, is 


“ “ . TUttX . 71tV 

w= X 1 X a mn sin sill 

m = 1 n = l & v 


(E4) 


where the origin is taken at a corner of the plate. Substi- 
tuting in equations (11) and (12) 

Jm(x )= sin — 


and 


t X . niry 
g n (y)= sin -f- 


Q—Qi 


and performing the integration over the whole plate (limits 
x=0, a; y= 0, b) gives the set of equations 


[Y , , 9 a 2 \ 4 , 12Z 2 mV 

a mn ^m 2 +n 2 p J +~^i 


lc x m 2 p (m 2 +n 2 p) -k v n 2 p (m 2 +n 2 p) J 


™ 2 +^ 2 p) 1 + 


321, ^ 

7T 2 6 3 (m 2 —p 2 ) (n 2 —q 2 ) 


(E5) 


where m=l, 2, 3, ... , w=l, 2, 3, ... , and p and q take 
only those values for which m±p and n±q are odd numbers. 

Equation (E5) represents an infinite set of homogeneous 
linear equations involving the unknown deflection coeffi- 
cients da. In order for the deflection coefficients to have 
values other than zero, that is, in order for the panel to 
buckle, the determinant of the coefficients of the unknown 
deflection coefficients must vanish. This determinant can 
be factored into two subdeterminants, one involving the 
unknown deflection coefficients a for which i±j is odd 
and the other involving those coefficients for which i±j is 
even. Buckling occurs, therefore, when either of the two 
subdeterminants vanishes. Only the buckling criterion in- 
volving the even subdeterminant is treated here. This 
criterion is 
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where 


M m 


= 32 Zi* [( /w2+ri2 ID' + ' ^ P ( m2+?i2 p) ~ Kn2 P ( m2+n2 p) ] 

Division of each column of the determinant in equation (EG) by the proper 

(**•»■ 

gives the simplified equation 

a n + ^13^1 + 9 ^ a 22 ^4+4 ^ a 3I ^9 + ^ #33^9 + 9 ^ . . . 


A 7 n 

0 

4 

9 


0 

A 7 13 

4 

5 


4 
9 

_4 

5 

AA 


3G 


_4 

5 

A 31 


0 

3G 

25 
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A 7 33 


= 0 


where 

T 2 b 3 

A mn 32 Avr 


(m 2 +n 2 
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7T 4 6 4 


1 2 Zrm*a* 

( m2+n2 p) 


7 •>« 7 

2 y—k x m- b2 —k y n- 


b A 


The vanishing of this determinant is the criterion for the 
symmetrical buckling of the shell. The same buckling cri- 
terion results from the use of the modified equation, as is 
shown in the following section. 

SOLUTION OF MODIFIED EQUATION 

The modified equation expressing the equilibrium of a 
curved panel under constant median-surface stresses in 


(EG) 


(E7) 


general form is 

. d 2 w . _ . d 2 w . d 2 w „ . 
DV W+ r 2 d/' + aJ bx> +2rt dxdy +avt dP = 0 (E8) 

Division of equation (E8) by T) and simplification of the 
result gives the following equation: 


12Z 2 4 d 4 w , 7 rd 2 w 


t 2 d 2 w 




Equation (E9) can be represented by 

Q 2 (w)= 0 
where Qo is defined as the operator 

i2Z* a> x 2 a 2 x 2 a 2 * 2 a 2 

v + p v ax: 4 p ax 2 r -■ * b 2 ax a ,y Kv p a p 


(E9) 

(E10) 
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By use of the Galerkin method and by use of the expression 
for w given in equation (E4), the following set of equations 
analogous to equations (E5) are obtained 


(m*+7^y+ - 


1 2 Z 2 m 4 o 


22k s w' 


7r'i 4 ( m 2 +?i- 2 p) 


m n /ii/ 


'(■ 


a ->\2 — k x m 2 ^ 2 —k,jTi- 


b* 


+ 


]C >- 


x 2 6 3 pi pi vt (m 2 — p 2 )(n 2 — <z 2 ) ' 


(Ell) 


where m=l, 2, 3, . . . , n=l, 2, 3, ... , and j) and q take 
only those values such that m±p and n±q are odd numbers. 

As in the case of the solution of Donnell’s equation, the 
stability determinant representing equations (Ell) can be 
factored into an even and an odd subdeterminant. The 
even one is 
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013 

022 

031 

033 

111 = 1 , 

n— 1 

N n 

0 

4 

9 

0 

0 

m= 1, 

ti=3 

0 

N tt 

4 
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0 

0 

m = 2, 

71=2 

4 

9 

4 

5 

N l2 

4 

5 

36 
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m= 3, 

n=l 

0 

0 

_4 

5 

iv 31 

0 

m = 3, 

71 = 3 

0 

0 

3G 

25 

0 

^33 


= 0(E12) 


The stability determinant (equation (E12)) obtained from 
the modified equation is identical with the simplified stability 
determinant (equation (E7)) obtained by use of Donnell’s 
equation. This identity holds for the odd as well as the even 
determinants. 

Although the stability determinants obtained by use of 
the two equations are identical and yield identical buckling 
loads, the determinant in equation (E7) consists of the 
/ . . a 2 V 

coefficients of a A i 2J rj 2 p) > whereas the determinant in 

equation (El2) consists of the coefficients of Accord- 
ingly, although the buckling loads found by the two methods 
are the same, the buckle patterns are different. Of the two 
buckle patterns the one found by the use of the modified 
equation is believed to be correct. This conclusion has been 
verified for the limiting case of a flat plate (Z= 0). 
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Positive directions of axes and angles (forces and moments) are shown by arrows 




Axis 


Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 


Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling _ . _ 

L 

Y >Z 

Roll. 


u 

V 

Lateral .. . 

Y 

Y 

Pitching 

M 

z — 

Pitch 

e 

0 

Q 

Normal. 

Z 

Z 

Yawing. 

N 

X >Y 

Yaw 


w 

r 


Absolute coefficients of moment 

a C -N_ 

Vm ~qcS '~ >n qbS 

(rolling) (pitching) (yawing) 


e - & 
Cl ~qbS 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.} 


4. PROPELLER SYMBOLS 


D Diameter 

p Geoinetric pitch 

p/D Pitch ratio 

V' Inflow velocity 

V 8 Slipstream velocity 

Thrust, absolute coefficient C T = 


T 

Q 


Torque, absolute coefficient C Q = 



Cg 


V 

T 

n 

pv?D* 


Q 

$ 

ptfD* 



P Power, absolute coefficient 

Speed-power coefficient = ^p !~2 
Efficiency 

Revolutions per second, rps 
Effective helix angle ^tan -1 ^^^ 


5. NUMERICAL RELATIONS 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 
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